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Abstract 

The cosmological remnants of a first-order phase transition generally depend on 
the perturbations that the walls of expanding bubbles originate in the plasma. Sev- 
eral of the formation mechanisms occur when bubbles collide and lose their spherical 
symmetry. However, spherical bubbles are often considered in the literature, in par- 
ticular for the calculation of gravitational waves. We study the steady state motion 
of bubble walls for different bubble symmetries. Using the bag equation of state, we 
discuss the propagation of phase transition fronts as detonations and subsonic or 
supersonic deflagrations. We consider the cases of spherical, cylindrical and planar 
walls, and compare the energy transferred to bulk motions of the relativistic fluid. 
We find that the different wall geometries give similar perturbations of the plasma. 
For the case of planar walls, we obtain analytical expressions for the kinetic energy 
in the bulk motions. As an application, we discuss the generation of gravitational 
waves. 



1 Introduction 

Cosmological phase transitions generically produce cosmic relics such as topological de- 
fects pQ , magnetic fields [2] , baryon number asymmetries [3] , inhomogeneities [U [5] , and 
gravitational waves [HI [7] . The walls of expanding bubbles usually play a relevant role in 
the mechanisms which generate these relics. Indeed, most of them depend on the distur- 
bance produced by the motion of bubble walls in the surrounding plasma. Moreover, the 
wall velocity itself depends on such hydrodynamics PEJlinillllinillSlIIlIISlIISlinias 
well as on the friction with the plasma [18]. The friction is determined by microphysics, 
i.e., the interactions of particles with the wall. The hydrodynamics is determined by the 
relativistic fluid equations. The latent heat that is injected at the phase transition fronts 
spreads out, causing reheating and bulk motions of the plasma. There are essentially 
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three propagation modes for the phase transition front. A detonation, which is super- 
sonic and is followed by a rarefaction wave, a subsonic deflagration, which is preceded by 
a shock front, and a supersonic deflagration, which is preceded by a shock and followed 
by a rarefaction wave. 

Several of the mechanisms which generate the aforementioned cosmological relics op- 
erate when bubbles meet. For instance, in order to produce gravity waves (GWs), bubbles 
must collide and lose their spherical symmetry, since a spherical source cannot generate 
gravitational radiation. The hydrodynamics of colliding bubble walls is complicated. The 
fluid velocity and temperature profiles during bubble collisions were studied using numer- 
ical simulations in Ref. [12]. In Ref. [13], a few configurations for the collision of shock 
fronts were considered analytically. Both studies were performed in 1+1 dimensions. In 
applications, further simplifications are needed. For instance, in the case of gravitational 
waves, the relevant quantity is the energy that is injected into bulk motions of the plasma. 
This is parameterized by the efficiency factor n, which is defined as the ratio of the kinetic 
energy in bulk motions to the released vacuum energy. In calculations of k, the interac- 
tion between bubbles is neglected, and spherical bubbles are assumed. This treatment 
implies two hypothesis, namely, that the motion of a bubble wall is not affected by the 
perturbations other bubbles caused in the plasma, and that the deformation of the bubble 
wall does not affect the energy transfer from the wall to the plasma. The first assumption 
should be correct for supersonic walls, which either do not have shock fronts preceding 
them or the shock fronts are very close to the walls. For subsonic walls, the main influence 
from other bubbles is the reheating of the plasma. This effect must be taken into account 
in a complete calculation of the phase transition [5j [191 120] • Regarding the dependence 
on the wall geometry, some progress can be made by comparing different bubble symme- 
tries. Of course, after a few collisions bubbles will take arbitrary forms. Nevertheless, 
considering a few specific symmetries will be useful to tell whether the disturbance of the 
fluid has a strong dependence on the wall geometry or not. Furthermore, it is particularly 
important to study the case of planar walls, for which analytical results can be obtained. 

The growth of plane, cylindrical, and spherical bubbles (equivalently, bubbles in 1, 
2, and 3 spatial dimensions, respectively) was considered in Ref. [TO] for the case of 
deflagrations. The fluid velocity profiles are quite different. In particular, for planar walls 
the shock wave preceding the phase transition front has a constant fluid velocity, whereas 
in the cylindrical and spherical cases the fluid velocity falls quickly in front of the wall. 
This is because in higher dimensions there is more room for the shock wave to carry away 
the injected energy. Nevertheless, the total amount of kinetic energy of the fluid must 
be a fraction of the released latent heat (another fraction goes into thermal energy), and 
there is no reason for this fraction (and for the ratio n) to have a strong dependence on the 
bubble wall geometry. For planar walls, the velocity profile of the fluid can be obtained 
analytically (see, e.g., [TO]), and one expects to find analytical formulas for quantities such 
as the efficiency factor as well. An analytical approximation for k was obtained in Ref. 
[2T] for small wall velocities. Recently, the efficiency factor was calculated numerically for 
spherical bubbles [22] in the whole wall velocity range. 

In this paper we address the issue of the effect of the wall geometry on the distur- 
bance caused by the motion of the walls. Thus, we study the hydrodynamics of spherical, 
cylindrical, and planar bubble walls. In particular, we show that the efficiency factor does 
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not differ significantly between the different wall geometries. For planar walls, we obtain 
the efficiency factor analytically for any wall velocity. We also discuss the generation of 
GWs in the electroweak phase transition. The plan of the paper is the following. In the 
next section we review the dynamics of the first-order phase transition. In section [3] we 
calculate the fluid profiles and the efficiency factor for the three bubble symmetries. In 
section H] we study in detail the case of planar walls. We discuss the different hydrody- 
namical modes and we derive analytical formulas for k. In section [5] we apply our results 
to the estimation of the gravitational wave signal from the electroweak phase transition. 
The conclusions are in section El 

2 Phase transition dynamics 

A cosmological phase transition occurs when the free energy of the model (i.e., the finite- 
temperature effective potential) depends on an order parameter (e.g., a Higgs field), such 
that the free energy density T ((f), T) has a minimum + (T) at high temperatures, and 
a different minimum 0_ (T) at low temperatures (we shall use a "+" index for variables 
in the high-temperature phase and a "— " index for the low-temperature phase). For 
a first-order phase transition, there is a temperature range in which the two minima 
coexist and a barrier in the free energy separates them. The critical temperature T c 
is that at which the two minima have the same free energy. Below T c , bubbles of the 
stable phase nucleate and grow. The nucleated bubble is a configuration cf) = <ft (r, t) with 
spherical symmetry, such that at the center of the bubble the system is in the low-T 
phase, whereas far from this point the system is in the high-T phase. Hence, we have 
(ft (r = 0) = 4>- (T) and 4> (r = oo) = (p + (T). There is a region, the "bubble wall", in 
which varies continuously from 0_ to <f>+. In this work, we shall assume for simplicity 
an infinitely thin wall separating the two phases. 

We are interested in the energy in bulk motions of the fluid which is caused by the 
moving walls. We assume that the plasma is a perfect relativistic fluid with four-velocity 
field m m = (7,7V), with 7 = l/yl — v 2 . The energy- momentum tensor is of the form 

T M " = (e+p) u^u" — pg^ u , (1) 

where e and p are the energy density and pressure in the proper system of the fluid element 
[23] . The energy density is given by T 00 = wy 2 —p= (e + pv 2 ) 7 2 , and the kinetic energy 
density is defined as ekm = T 00 (v) — T 00 (0). Therefore, we have ekin = T 00 — e = wv 2 ^ 2 . 

All the thermodynamical quantities can be derived from the free energy densities 
T+ (T) = T (0 + (T) , T) and T- (T) = T (</>_ (T) , T) for each phase. Thus, the pressure 
is given by p = —J 7 , the entropy density by s = dp/dT, the energy density by e = Ts — p, 
and the enthalpy by w = e + p = Ts. At the critical temperature, the pressure in the two 
phases is the same, i.e., p + (T c ) = p_ (T c ). However, other quantities such as the energy, 
entropy, and enthalpy are different even at T = T c . The latent heat is defined as the 
energy density difference L = Ae (T c ) = Aw (T c ) = T c As (T c ). 

It is useful to consider a simplified model which exhibits the general features of a 
phase transition. This allows in particular to obtain results which depend on a few 
parameters. Then, the results can be applied to realistic phase transitions by calculating 
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these parameters in specific models. Therefore, we shall use the bag equation of state 
(EOS), 

e + = a + T 4 +e, e_ = a„T 4 , , . 

p+ = |a + T 4 — e p_ = |a_T 4 , 

with e > and a + > a_ > 0. This system has two components, namely, a vacuum energy 
density e and a radiation energy density aT 4 . The condition a + > a_ implies that some of 
the radiation degrees of freedom disappear after the phase transition. The vacuum energy 
density is positive in the high-temperature phase and vanishes in the low-temperature one. 
In this model, the critical temperature is given by T c = (35/Aa) 1 ^ 4 , where Aa = a + — a_, 
and the latent heat is related to the vacuum energy density by L = 4e. In both phases, 
the speed of sound is given by 

c 2 s = dp/de = 1/3. (3) 

We shall assume for simplicity a stationary state in which the wall moves with a 
constant!]] velocity v w . We aim to consider three kinds of bubble wall geometry, namely, 
a spherical, a cylindrical and a plane wall. In this approximation, at time t the wall is 
at a distance Rb = v w t from a point, axis, or plane, and the volume of the bubble is of 
the form V& = CjB? b + / (j + 1), where j = 2, 1, or for the spherical, cylindrical or planar 
case, respectively (the factor Cj will cancel out in our calculations). The kinetic energy 
in bulk motions of the fluid is given by E^ n = Cj J °° wv 2 ^ 2 R^ dR. The efficiency factor is 
defined as the ratio of the kinetic energy to the released vacuum energy, 

k = E kin / (eV b ) . (4) 
J + l 



Thus, we can write 



ea +1 Jo 



wv'Yt'dt, (5) 



where £ = R/t and £ w — Rb/t — v v 



3 Fluid profiles and kinetic energy 
3.1 The fluid equations 

The fluid equations are obtained from the conservation of the energy-momentum tensor 

If we denote by r the distance from the symmetry point, axis or plane, and t the time 
from nucleation, Eqs. (EJ) become [10] 

d t [{e + pv 2 )-i 2 ]+d r [(e + p)-{ 2 v\ = - 3 - [(e + p) -y 2 v] , 

d t [(e + p) 1 2 v]+d r [{ev 2 +p) 1 2 ] = - J -[(e + p)^v 2 ]. (7) 

a real phase transition, even after reaching the stationary state, the wall velocity may vary due to 
the adiabatic cooling of the universe and the release of latent heat (see, e.g., [5l H"9l [24] ) . 
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Since there is no characteristic distance scale in the problem, it is usual to assume the 
similarity condition, namely, that e,p and v depend only on the combination £ = r/t. 
Thus, we have 



(£-«)- = j- + 7 2 (i-<K 

w 4 
w 



(8) 



where a prime indicates derivative with respect to £. The pressure and energy density are 
further related by the equation of state. According to Eq. (j3J), we have p' = c 2 s e' . Using 
this relation, Eqs. dHJ) can be combined to obtain the central equation for the velocity 
profile [SI [HUE]. We obtain 



v . 



where 



From Eqs. flS]) we also obtain the equation for the enthalpy profile, 

1 



w 
w 



+ 1 ] /ryV, 



(9) 



(10) 



which is readily integrated [7]. For c 2 s = 1/3 we have 



w b 
— = exp 



4 7 2 /i(e,t;)^ 



'12) 



Notice that, in this model, the equations for the velocity and enthalpy profiles are the 
same in both phases. 



3.2 Discontinuities 

3.2.1 The phase transition front 

In the £ axis, the bubble wall is at £ w = v w . The enthalpy and other quantities are 
discontinuous at and so will be the fluid velocity. In the reference frame of the wall, 
the fluid comes from the high-T phase with a velocity v + , and goes out into the low-T 
phase with a velocity t>_. The incoming and outgoing flow velocities are related by the 
conservation of T^ v across the wall 



Equivalently, 



v + v. 



Pi 



P- 
e_ 



v+ 

V- 



e- +P+ 
e++P-' 



(13) 
(14) 

(15) 
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Notice that these equations do not depend on j. This is because the surface of disconti- 
nuity is locally planar. In the reference frame of the bubble center, the fluid velocities on 
each side of the wall are given by v± = \i (£ w , \v±\). 

According to Eq. ([Hj) , v + and v_ have the same sign. Indeed, the fluid velocities 
in the system of the wall must be both negative. Using the bag EOS, Eqs. ffT5l) can be 
combined to obtain a relation between v + , t>_ and the parameter 

a + = r^i • ( 16 ) 

We can solve, e.g., for v + as a function of t>_ and a + [8], 



+ <t ) ± V ( ¥ + 6^1 j' + C 1 + «+) («+ " 1/3) 
1 + a + 

Two kinds of hydrodynamical processes may occur at the phase transition front, corre- 
sponding to the + and — signs in Eq. (I17p . namely, a detonation, for which the incoming 
flow is supersonic (\v+\ > c s ) and faster than the outgoing flow (|i>_| < |f+|), and a 
deflagration, with \v + \ < c s and > In either case, the incoming velocity \v + \ 

has an extremum at \v-\ = c s , namely, a minimum for detonations and a maximum for 
deflagrations. A process with |i>_| = c s is called a Jouguet detonation or deflagration. In 
this case we have \v + \ — Vj (a+), with 



< («+) = TfTT- v • (18) 

V3 (1 + a+) 

Hence, for detonations we have c s < v j ct (a+) < \v + \, and for deflagrations \v + \ < 
Vj ei (a + ) < c s . The hydrodynamical process is called weak if the velocities v + and V- 
are either both supersonic, or both subsonic. In such a case, |t>_| lies between \v + \ and 
c s . Otherwise, one of the two velocities is subsonic and the other one is supersonic. Such 
hydrodynamical process is called strong. 

3.2.2 Shock fronts 

Discontinuities in the same phase may also be needed to satisfy the boundary conditions. 
Such discontinuities are called shock fronts. We shall use the index 1 for fluid variables 
behind the shock front at £ = £ s h, and the index 2 for variables in front of the shock. The 
EOS is the same on both sides of the discontinuity, and Eqs. (IT5|) trivially give 

1 Vl 3T 2 4 + Tf 

^ 2 = 3' - 2 = *TTTTi> (19) 

where v i and v 2 are the (negative) fluid velocities in the shock frame. 

When dealing with discontinuities, entropy considerations can be useful for discarding 
possible processes (see, e.g., [9]). Consider a portion of the fluid which passes through 
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a discontinuity surface. Requiring the entropy of the fluid to increase, one obtains the 
condition Sit>i7i > s 2 v 2 ^ 2 . Using s = w/T and Eq. ffT4"l) . this condition becomes 

T 2 /T 1 > 7l / 72 . (20) 

For the bag EOS, we have W2/1V1 = (T 2 /Ti) A , and we may insert the inequality f[2T)j) back 
in Eq. ([H]) to obtainl 

v 2 {l-v 2 2 )<v l {l-v 2 1 ). (21) 
Using the relation V\v 2 = 1/3, this condition becomes (y\ — 1/3) 3 < 0, which implies 

\vx\ < 4= < H- (22) 

In the frame of the bubble center, the fluid velocities on each side of the shock front 
are given by = fi (£ s h, l^i^i)- According to the condition f[2"2"j) . we have 

V! > v 2 . (23) 

Hence, the fluid velocity must have a negative jump. As a consequence, V\ cannot vanish 
(otherwise, we would have v 2 < 0, but, in this reference frame, the fluid velocities are 
positive or zero). On the other hand, we may have v 2 = 0. In this case, the velocity of 
the shock front is given by £ sh = — v 2 , and the first of Eqs. f|T9l) gives 

vi = (24) 

Equivalently, 



v 1 / / vi \ 2 , 1 



which implies that the shock is supersonic. 



3.3 Kinds of solutions 

Equation (jH]) can be solved numerically for the spherical and cylindrical cases, and analyt- 
ically for the planar case. General solution curves for the three symmetries can be found 
in Ref. [10J (for the planar case, see section H] below). The velocity profile of the fluid 
must fulfil the discontinuity conditions at the bubble wall, and is constructed by matching 
different solutions of Eq. The boundary conditions are that the fluid is at rest far 
ahead of the phase transition surface (where no information of the phase transition has 
arrived yet) and far behind that surface (near the center of the bubble). Far in front of 
the phase transition front, the temperature is still Tjv- Therefore, the boundary condition 
for the enthalpy density is that it takes the value 

4 

w N = -a + T^ (26) 
2 Similarly, for the phase transition front we obtain u+(l — vi) < (a_/a + )v_ (1 — v 2 ^). 
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far in front of the bubble wall. 

Not all of the aforementioned hydrodynamical processes will be realized in a phase 
transition. It is known that strong detonations are not possible, since they cannot satisfy 
the boundary conditions [231 El Q3] . It has been argued in Ref . [5] that weak detonations 
are also impossible, like in the case of chemical burning [23]. As a consequence, we would 
only have Jouguet detonations, and the velocity v w = Vj would be completely determined 
by hydrodynamics and would not depend on microphysics. However, it has been shown 
[2] that this is not true in the case of phase transitions, where the situation is similar 
to that of condensation discontinuities [23] rather than chemical burning. Consequently, 
only strong detonations are forbidden. 

Strong deflagrations do not seem to be realized either. In Ref. [H] it was argued that, 
like in the case of chemical burning, they are forbidden by entropy considerations [23] . 
However, in Ref. [13], it was shown that this proof is not valid for cosmological phase 
transitions. Nevertheless, a mechanical instability argument against strong deflagrations 
[23] seems to be valid also for cosmological phase transitions [15] . Numerical calculations 
[T2"l [15] support this assertion. On the other hand, in Ref. [15] supersonic Jouguet 
deflagrations were shown to exist. 

As a consequence, three kinds of solutions seem to be realized in nature, namely, 
weak detonations, subsonic weak deflagrations, and supersonic Jouguet deflagrations. In 
section HI for the planar case, we argue that these are the only possible solutions (the 
argument is simpler for the planar case, but can be straightforwardly generalized to the 
other symmetries). In the rest of the present section we describe these solutions and cal- 
culate their profiles and the efficiency factor (the spherical case has been studied recently 
in Ref. [22]). 

3.3.1 Detonations 

For detonations we have \v + \ > v j et , and the wall moves supersonically with respect to 
the fluid in front of it. Thus, outside a detonation bubble the fluid has not yet been 
perturbed. Hence the fluid velocity vanishes and the temperature is still that at which 
the bubble nucleated, i.e., v + = and T + = Tn. Therefore, we have 

v + = -£ w , a + = a N , (27) 

where 

a N = -4^- (28) 

As a consequence, we have £ w > Vj et . The velocity profile of a detonation is shown in 
Fig. [1] (left) for the three wall geometries. The bubble wall is followed by a rarefaction 
wave which ends at £ = c s . The velocity profile is determined by the boundary condition 
v{£w)=v-, with 

iL = fi (&,>_|). (29) 
The velocity t>_ is given by the inverse of relation (fTTI) . 

m - + ) ± \Z(^^ + ) 2 -I <*> 
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Weak detonations (i.e., |i>_| > c s ) correspond to the + sign in Eq (1301) and strong deto- 
nations (\v-\ < c s ) correspond to the — signal. The two branches match at the Jouguet 
point \v + \ = Vj et , \v-\ = c s . As we mentioned, the process must be a weak (or, at most, 
Jouguet) detonation for compatibility with the boundary conditions. 

The enthalpy profile is given by Eq. (1121) with the condition w (C, w ) = The value 
W- just behind the wall is related to w + = wjy through Eq. (TH|) . 



w. 



I I 2 



-W N . 



(31) 



2^,2 



The efficiency factor is obtained by integrating the kinetic energy density e^m = wV'j 
Figure [T] (right) shows the kinetic energy density profile for the three wall geometries 




Figure 1: Left: the fluid velocity profile of a detonation with v w = 0.8 and o-y = 0.1 for a 
spherical wall (solid), a cylindrical wall (dashed), and a planar wall (dotted). Right: the 
corresponding kinetic energy density profiles. 



3.3.2 Subsonic deflagrations 

For £ w < c s , the wall is preceded by a shock front at £ s h > c s . Behind the wall the fluid 
velocity vanishes. Hence, we have V- — and (since v + > 0) \v + \ < \v-\ — £ w < c s , i.e., 
the hydrodynamical process is a weak deflagration (in the limiting case £ w — c s , we have 
a Jouguet deflagration). The fluid velocity vanishes also beyond the shock front. Thus, 
we have v 2 = and a 2 = a^- 

The velocity profile of the shock wave between £ w and £ s h is given by a solution of Eq. 
dSJ) (see Fig. [2]). One can choose as boundary condition the value of the velocity v + in 
front of the wall, or the value v\ behind the shock. The velocity v + is given by \i (£ w , \v+\), 
where v + is given by Eq. (TTTj) as a function of a + and of f_ = —v w . However, in the 
present case a + does not equal a^, since the shock wave reheats the fluid in front of the 
wall. The temperatures at each end of the shock wave are related by 

a + = a 1; {62) 

w + 

3 In Eq. (f5U|) . |i7 J is real only for \v+\ > Vj ct (corresponding to detonations) or \v + \ < v d j cf (corre- 
sponding to deflagrations). 
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with w\/w + given by Eq. (|T2|) . 



Wi 

— = exp 

w + 



vi 



v+ 



(33) 



The fluid velocities v + and 5i are related by the fluid equation (j^J). Let us denote t> (£; v\) 
the solution with boundary condition v (£ s h) = t>i, with £ sh depending on V\ through Eq. 
(I25|) . Then, evaluating at £ = we obtain 5 + as a function of 5i and £ w , 

v+ = v(£ w ;vi) . (34) 

Equations (I321I34P give a + as a function of a%, £ w and v%. The variables ai and a.i = a at 
on each side of the shock discontinuity are related by Eqs. ( Tl9|) with t> 2 = — £ s h- We 
obtain 

«i = Q V 2 _t ( 35 ) 
Equivalently, using Eq. (|25|) . we have 



ai = ^ ^3 + 5^ 2 - 401^3 + a*. (36) 

Thus, we have o;_)_ clS cl function of «at, and either v\ or £ s h. On the other hand, the 
variables at the wall discontinuity are related by Eq. fTTTj) or Eq. fl30|) or, equivalently, by 

"+ = 7+ ^+ + ^ - v+v- - • (37) 

In the present case, we have t>_ = — £ w , v + = (v + — £ w ) / (1 — v+£ w ), and we obtain 

~2 ~ 

a + = ^± (2^*5+ + 1 - 3i£) . (38) 

Inserting Eq. ( 158]) in Eq. f )32|) . we eliminate a + and we can solve for v\ as a function of £ w 
and «Ar. Searching for v\ numerically, implies the evaluation of all the quantities which 
appear in Eqs. (|32H38p . for several values of V\. Such evaluation involves numerically 
solving the differential equation which gives v + as a function of v\ and then performing 
numerically the integral in Eq. (133]) . As we shall see in section HJ in the planar case Eqs. 
(I32H38P reduce to a single, algebraic equation, which can be solved analytically. 

Once the value of v% is found, one can compute the velocity and enthalpy profiles and 
perform the integral of the kinetic energy density to obtain the efficiency factor. The 
enthalpy profile is given by Eq. (fl2"j) and determined by the condition w (£ s h) = Wi, where 
w\ = (at/v/cKi) w n, with a^/ax given by Eq. fl35|) . Thus, we have 



*" = W^f' (39) 
Figure [2] shows the velocity and kinetic energy density profiles for a subsonic deflagration. 
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Figure 2: Left: the fluid velocity profile of a deflagration with = 0.4 and = 0.1 for 
a spherical wall (solid), a cylindrical wall (dashed), and a planar wall (dotted). Right: 
the corresponding kinetic energy density profiles. 



3.3.3 Supersonic deflagrations 

For c s < £ w < Vj ct , the phase transition front is preceded by a shock front at £ s h > £ w , and 
is followed by a rarefaction wave solution which vanishes at £ — c s . Therefore, both ve- 
locities v + and u_ are non- vanishing (see Fig. [3]). The hydrodynamic process is a Jouguet 
deflagration, i.e., \v + \ < |t>_| = c s [15J. The wall velocity £ TO = (u_ + c s ) / (1 + iLc a ) is 
supersonic and depends on the value of tL. For = this solution matches the "ordi- 
nary" deflagration considered before. As £ w approaches v j et , the shock wave gets shorter, 
i.e., £ s h —> £ w , and the profile matches that of the detonation considered before. 

The calculation of the boundary value v± for the shock wave profile is very similar 
to that of the ordinary deflagration. Indeed, Eqs. (I32H37P hold in this case. Since the 
deflagration is now Jouguet, we have the condition t>_ = — c s instead of i>_ = — v w . As a 
consequence, Eq. (j3"Tj) becomes a + = 7+ (v + + , and Eq. fl38|) gets replaced by 



a + = — - — 



(l - V3v w - v+(v w - V3)Y . (40) 



From Eqs. (1521 and ([40]) we eliminate a + and we obtain Vi as a function of £ w and 
as before. The enthalpy profile of the shock wave is determined by the value w\ at the 
shock, given by Eq. (139|) . 

The profile of the rarefaction wave is similar to that of the detonation, with the 
boundary condition v = 5_, with 5_ now given by 

& ~ C s (41) 



f Sw^s 

The enthalpy profile behind the wall is determined by the condition w (£ w ) = The 
value of the enthalpy just behind the wall is now given by 

w. = (2/v / 3)|^ + |7>+, (42) 
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with v + given by 

v+ = {v + - U / (1 - S+U , (43) 

and v + and w + are obtained from the shock wave profile. Figure [3] shows the profile of a 
supersonic deflagration for the three geometries. 
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Figure 3: Left: the fluid velocity profile of a deflagration with v w = 0.7 and = 0.1 for 
a spherical wall (solid), a cylindrical wall (dashed), and a planar wall (dotted). Right: 
the corresponding kinetic energy density profiles. 



3.4 Energy injected into the plasma 

The latent heat released at the phase transition fronts spreads out in the plasma. Part 
of this energy causes reheating, and another part causes bulk motions. Some of the 
consequences of the phase transition will depend on the thickness of the plasma shell 
where the energy is concentrated. The issue of reheating was addressed, e.g., in Refs. 
[21 [H] for the electroweak phase transition. The size of the regions of reheated plasma 
affects the dynamics of the phase transition and, as a consequence, the baryogenesis 
mechanism. Here we shall focus on the energy in bulk motions of the plasma. 

The bubble walls set the fluid moving forward with velocities v + and In the left 
panels of Figs. [U [2J and [3] we observe that these velocities are very similar for different 
wall geometries. However, away from the wall the fluid velocity decays faster for more 
symmetric bubbles, since in that case there is more room for the energy to get distributed. 
For illustrative purposes, we define the shell where kinetic energy is concentrated as the 
region around the wall where the kinetic energy density remains higher than half the 
maximum on each side of the wall. Thus, the thickness is given by 5£ = £ + — where 
£ + and are determined by the condition ey m (£+) = 0.5w±5±7±. We plot the value of 
5£ in Fig. |U The shape of the curves is different from those given in Ref. [22] for the 
spherical case, since we have defined <5£ differently. Nevertheless, the general structure 
is similar: the thickness is larger for subsonic deflagrations than for detonations, and is 
quite small for Jouguet solutions. We also see that the energy is more widely distributed 
in the planar case, especially for weak deflagrations. 

Some of the cosmological remnants of the phase transition do not depend on the 
thickness of the region of perturbed fluid around the wall. This is the case, e.g., when the 
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Figure 4: The thickness <5£ of the region around the wall inside which the kinetic energy 
density decreases to a half of its maximum value, for oln = 0.1 as a function of v w . Solid 
lines correspond to spherical bubbles, dashed lines correspond to cylindrical bubbles, 
and dotted lines correspond to planar walls. Subsonic deflagrations are plotted in blue, 
supersonic deflagrations are in black, and detonations in red. 

generating mechanism is based on the turbulence produced by the colliding bubble walls. 
Eddies are formed at all size scales up to the size of the largest bubbles. 

It is usually assumed that detonation walls cause a stronger disturbance of the fluid 
than deflagration walls, since detonations have higher velocities. For weak detonations, 
however, as pointed out in Ref. [25], the higher the wall velocity, the smaller the distur- 
bance of the fluid. Thus, the strongest disturbance is caused by the Jouguet detonation 
(which is the case usually considered in the GW literature). However, this is just a lim- 
iting case; real detonations are generally weak. Furthermore, as noticed in Ref. [2T] . 
ordinary deflagrations may cause important perturbations in the fluid if they are close to 
the Jouguet limit v w = c s . In Fig. [5] we give the efficiency factor k as a function of v w for 
several values of a at. We calculated k numerically for the spherical and cylindrical cases, 
and analytically for the planar case (see section H]). For the spherical case (solid lines), 
our results agree with the numerical fits provided in Ref. [22J. As expected, for fixed a/v, 
the efficiency is larger for stronger solutions (i.e., solutions which are closer to Jouguet 
processes). Thus, for subsonic deflagrations the efficiency factor increases with the wall 
velocity, whereas for detonations k decreases with v w . The efficiency peaks for supersonic 
deflagrations, which are Jouguet processes. It is interesting to notice that even subsonic 
deflagrations can give larger efficiency factors than detonations. 

We see that the different geometries give in general similar values of k, except for 
the case of small wall velocities (v w < 0.2). This case will not be interesting in general 
for GW generation, since the efficiency factor is small. For faster walls, the three curves 
are quite close; they only separate somewhat for supersonic deflagrations (in Fig. [51 the 
differences at the top of the curves are less than a 20%). In any case, it is clear from Fig. 
Othat the difference between geometries will always be at most an O (1) factor, except 
for uninterestingly small values of k. 
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Figure 5: The efficiency factor k as a function of v w for a N = 0.01,0.03,0.1,0.3. Solid 
lines correspond to the spherical case, dashed lines correspond to the cylindrical case, and 
dotted lines correspond to the planar case. Subsonic deflagrations are plotted in blue, 
supersonic deflagrations are in black, and detonations in red. 

4 Planar walls: analytic results 

In this section we calculate analytically the kinetic energy density and the efficiency 
coefficient for planar walls. In this case the "bubble" consists of two planar walls at 
positions x = v w t and x = —v w t, and is equivalent to a bubble in 1+1 dimensions. The 
system is symmetric under reflection through a plane, and we need only consider the 
wall moving to the right. The absence of a length scale in the problem implies that the 
profiles depend only on £ = x/t. The solutions of the fluid velocity equation for the 
planar case are well known. For j = 0, Eq. ([9]) implies either that = or that 

/i (£, i>) = ±c s . The latter implies that v — \i (£, ±c s ). These solutions are shown in Fig. 
[H]for £ > and v > 0. We need only consider that quadrant since the wall at £ w = v w 
sets the fluid moving forward (the reflected profiles around the opposite walfl at £ = —v w 
are constructed with the solutions for £ < and v < 0). As we shall see, the solution 
v — fi (£, — c s ) will not take part in the fluid profile. Indeed, the possible values of the 
fluid velocity will be those below the curve v = £ (in dots in Fig. [6j) due to the fact that, 
at the wall, the fluid velocity fulfils v± < £ w (since in the wall frame we have v± < 0). As 
a consequence, the physical solutions are either the constants or the "rarefaction" 



The solution t> rar is positive only for £ > c s . Between £ = c s and £ = 1, f ra r(0 grows 
monotonically from v = to v = 1. 

4 More complicated profiles will arise if one considers two bubbles nucleated at a distance d [13]. In 
particular, this separation introduces a new length scale in the problem. 



^rar (£) 



i - c s 

l-c s £' 



(44) 
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Figure 6: Solutions to the fluid velocity equation for the planar case (solid lines). The 
dotted line is the curve v — £. Physical solutions are below this curve. The dashed line 
corresponds to the value of the fluid velocity v\ at the shock front. 

The enthalpy is readily obtained as well. For the case v' = the enthalpy is a constant. 
For the case \i = c s , the integral in Eq. ( 1T2|) is simple; given the condition w(£ a ) = w a , 
we obtain 



The velocity profile must be constructed by matching solutions in such a way that 
the discontinuity conditions at the bubble wall and the boundary conditions are satisfied 
(shock discontinuities may also be needed). Before obtaining the analytical results, we 
shall examine all the possible profiles and argue that only the three kinds of solutions 
considered in section [3] are acceptable. The argument is simpler for the planar case, but 
the generalization to the other wall geometries is straightforward. 

4.1 Hydrodynamic processes and fluid profiles 

The fluid velocity must vanish far in front of the wall, where no signal from the phase 
transition front has arrived yet. Since there are no decreasing solutions, we see that, no 
matter what the position £ w of the wall is, the fluid velocity must have a jump from v > 
to v — 0. This discontinuity can be either at £ = £ w or at some point £ s h between £ w and 
1. 

In the former case, we have v — for all In particular, v + = and, therefore, 

£ w = \v + \. Besides, the fluid velocity behind the wall must be nonvanishing, i.e., tL > 0. 
Thus, we have tL > v + and, hence, \v + \ > |t>_|. Therefore, the hydrodynamical process 
is a detonation. As a consequence, the wall is supersonic, since £ w = \v+\ > Vj et > c s . 
The fluid velocity must vanish also (by symmetry) at the bubble center. According to 
Eq. (I23p . the velocity cannot have positive jumps in the same phase. Therefore, 
must grow continuously from to £L. The only possibility is that v = for £ < c s and 
v = frar (£) for £ > c s . At some point between c s and £ w , the solution t> rar (£) may be 
matched again to a constant v = or it may continue growing until £ = £ w . In any case, 




(45) 
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we have £L < v TaT (£ w ) = fj,(£ w ,c s ), which implies that > c s . Hence, the detonation 
can only be weak or Jouguet. 

If the fluid velocity does not vanish in front of the wall (i.e., v + > 0), then it is clear 
from Fig. E]that the profile must have a shock discontinuity at some point £ s h > £ w , so 
that v = for £ > £ s h. The fluid velocity v\ at the shock front is given by Eq. (|24j) 
and shown in a dashed line in Fig. [6j We see that the solution t> rar lays completely on 
the right of the curve of t>i(£ s h)- Therefore, between £ TO and £ s h the solution must be a 
constant, u = tL = V\. 

Behind the wall, v can grow continuously from only if the wall is supersonic. 

If the wall is subsonic, then the fluid velocity must be v = for £ < £ w (otherwise we 
would need a positive jump). In this case we have V- < v + , which implies \v + \ < \v-\, 
and the process is a deflagration. Furthermore, since L = we have = £ w < c s , 
i.e., the deflagration is weak or, at most, Jouguet. 

If the wall is supersonic, we still have solutions for which v = behind the wall. In 
this case, the condition u_ = implies a strong deflagration (|i>_| = t; w > c s ). However, 
numerical simulations indicate that strong deflagrations are unstable [T2| [T5]. Notice that, 
since £ w is now > c s , we may have a non- vanishing fluid velocity u_ behind the wall, as 
in the case in which v + = 0. Like in that case, we have the condition v_ < t> rar but, 
instead of v + = 0, we now have v + = V\ (£ s h) > v Tar (£ w ) (see Fig. E]). Hence, we have 
S_ < v + and, thus, |f + | < Therefore, the present case is again a deflagration, 

not a detonation. Now V- can take any value, with the only condition {L < v TSir (^ w ), 
which implies that |i>_| > c s , i.e., the deflagration must be strong or Jouguet. Of all these 
solutions, though, one expects that the stable one will be that which causes the smallest 
perturbation of the fluid, i.e., the Jouguet deflagration \v-\ = c s . This is supported by 
numerical simulations. As we shall see, the supersonic deflagration matches the detonation 
solution at £, w = Vj et . 

Thus, a subsonic phase transition front always propagates as a weak deflagration and 
is preceded by a shock wave. A supersonic phase transition front is always followed by 
a rarefaction wave, and may propagate either as a Jouguet deflagration preceded by a 
shock front, or as a detonation, without a shock wave. 



4.2 Detonations 

For detonations, the fluid velocity is given by v — t> rar (£) between £ = c s and a certain 
£o < an d by v = iL between £ and £ w (see Fig. [[]). The matching condition 
^rar (£o) = determines the value of £o as a function of £L, 

io = — • 46 

1 + f_C s 

The velocity tL is given by Eqs. (125)1 and ( 15U|) as a function of a + = and f + = — £ w . 
The enthalpy is a constant w = w_ for £ < £ < £ tu , where w_ is given by Eq. (T5TT) as a 
function of and Between c s and £o, the enthalpy is given by Eq. (Jl5|) . with the 
condition io(£o) = Inserting the velocity profile (jSJ), we obtain 

(\ r 1 • ry'^ 
w = w_ — . (47) 
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Using Eqs. ( 147]) and OS]), the efficiency coefficient (j5]) is given by 
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The change of variable x = (l + £)/(l — £) leads to the simpler expression 

2 



(46 



(49) 



fix. 



(50) 



This integral can be expressed in terms of the hypergeometric function 2 F\ [2E]- We 
obtain 

1= \ if (&) - / (cs)) , (51) 
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(52) 



4.3 Subsonic deflagrations 



The profile for subsonic deflagrations is very simple in the planar case. The fluid velocity 
is a constant v = v + = V\ between £ w and £ s h, and vanishes outside that region. Thus, 
we have w + = W\ and a+ = a\, and Eqs. ( 1521155]) give an algebraic equation for V\. The 
equation is simpler in terms of £ s h, 



(3d - 1) 2 + U (3d - 1) = §<WtS, 



(53) 



Solving for £ s h as a function of «tv and v w amounts to finding the roots of a quartic 
polynomial. The algebraic expressions for the solutions are quite cumbersome and we 
shall not write them down. Only one of the four solutions gives £ s h > c s . The integral in 
Eq. §5§ is trivial since v is a constant, and the efficiency factor is given by 



1 ^l~2~2/ c c \ 



(54) 



where v\ is given by Eq. (124]) as a function of £ s h, and w\/e is given by Eq. ( 139]) as a 
function of £ w and ajv- 



5 See Eqs. 3.194-1, 9.131-1 and 9.137-2 of Ref . 
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4.4 Supersonic deflagrations 



In this case the shock wave in front of the wall is similar to that of subsonic deflagrations. 
We have again v + — vi, w + = Wi, and a + = a\. From Eqs. (I3"5"j) . (|4"0"|) . and (l2"4"|) we obtain 



7„ 



£ sh (l - v 7 ^) 



9 Q! TV 

4^' 



(55) 



which trivially reduces to a quadratic equation. The solution is 
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where x = a/3 — + y/ajv/'jw The rarefaction wave behind the wall is given by the 
solution t> rar (£). In the Jouguet case, CL is given by Eq. (l4"Tj) . which implies that £ — £«)■ 
The efficiency factor is given by 
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(57) 



The value of w_ in the Jouguet deflagration case is given by Eq. (|42[) as a function of 
w + = Wi, and depends on v + = {v\ — £ w ) / (1 — with ui given by Eq. pijl . The 

value of tui/e is again given by Eq. (139]) as a function of £ w and a at. 



For & 



,det, 



ajy), Eq. (|56|) gives £ s h = £ w [this is more easily checked by setting 



£sh = £w in Eq. (155]) ]. Hence, at the Jouguet detonation velocity, the shock disappears 
and the profile for the supersonic deflagration matches the profile for the detonation. 



5 Gravitational waves from real detonations and de- 
flagrations 

In section [3] we have discussed the disturbance that phase transition fronts cause on the 
plasma. In the present section we consider a particular consequence of such disturbance, 
namely, the generation of gravity waves. As we have already mentioned, gravitational 
radiation can only be produced once bubbles collide and lose their spherical symmetry. 
In fact, the "bubble collisions" mechanism P [23 ES] is based on the envelope approxi- 
mation [B] , which consists of taking into account only the motion of the uncollided walls. 
The thickness of the shell in which the energy of the fluid is concentrated is relevant 
for this mechanism, and the envelope approximation assumes that the energy concentra- 
tions are infinitesimally thin. On the other hand, at a cosmological phase transition the 
Reynolds number is large enough for bubble collisions to cause the onset of turbulence [7] . 
Turbulence turns out to be a more effective source of gravitational radiation than bubble 
collisions [7J [29j [30] . In an electrically conducting fluid and in the presence of magnetic 
fields, turbulence develops in a completely different way. This gives a third mechanism 
for generation of GWs in a phase transition (see e.g. [301 [3Tj). 
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The energy density of gravitational radiation is usually expressed in terms of the 
quantity 

h flaw if) = — -r: 7> 58 

Pc dlogf 

where pew is the energy density of the GWs, / is the frequency, and p c is the critical energy 
density today, p c = 3Hq/8ttG, with H = 100 frkrns - Mpc -1 . The GW spectrum depends 
on the details of the phase transition and on the generating mechanism. Nevertheless, 
the peak frequency f p is generally determined by the typical length scale of the source. 
In a first-order phase transition, the latter is the bubble size d, which is proportional to 
the duration At of the phase transition, d ~ v w At. The time At is in turn a fraction of 
the Hubble time. Once redshifted to today, the peak frequency is roughly given by 

/„ ~ 10- 2 mHz ^ 1 T * (59) 
Jp At lOOGeV v ; 

where if* and T* are the Hubble rate and the temperature at the moment of the phase 
transition. The sensitivity peak of the space interferometer LISA is expected to be 
h 2 Q ~ 10~ 12 at a frequency / ~ lmHz. Quite interestingly, GWs produced at the 
temperature scale of the electroweak phase transition, T* ~ lOOGeV, will have a charac- 
teristic frequency around f p ~ lmHz for At ~ 10~ 2 if which is a possible value for the 
duration of the phase transition. This motivated the investigation of the GW signal from 
the electroweak phase transition [251 E21 [33] . 

For simplicity, we shall consider only GWs from turbulence, for which we shall use the 
analytic approximation obtained in Ref. [SU] (we have checked, using the fit given in Ref. 
[28] . that the intensity resulting from bubble collisions is an order of magnitude smaller). 
The GW energy density Q p at the peak frequency depends on the length scale d ~ v w At. 
Thus, we have 

9n R 2 f At\ 2 2 2 j 4 for v < 1/2, 
^ * 32^ UrV " ^ I I/- 2 for v > 1/2, ( 6 °) 

where Qr ~ 5 x 10~ 5 is the radiation energy density parameter = Pr/ p c today, and 
v is the characteristic eddy velocity, defined by v 2 = jna N . Equation (p}|) should be 
valid both for detonations and deflagrations, although the time At must be calculated 
differently in each case. 

Notice that fl p is proportional to and to k 2 . The wall velocity is determined by 
hydrodynamics and by microphysics (see, e.g., [12]). As a consequence, v w depends on 
a friction parameter 77 as well as on the nucleation temperature Tjy. Different approxi- 
mations have been used for the friction force (see, e.g., [13(22]). For generality, we shall 
leave the result expressed in terms of v w . We do not expect v w to depend significantly 
on the geometry. In the case of deflagrations, the dependence of v w on T/v is affected by 
the shock wave, which depends on the wall geometry. In the case of detonations, v w is 
completely determined by the discontinuity equations (fT5j) and the friction, and does not 
depend on the geometry at all. In any case, as explained before, there is no reason to 
assume any particular symmetry after bubbles collide, and we shall use the planar wall 
results. 
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In Fig. [7] we plot the peak amplitude fl p of GWs from turbulence at the electroweak 
phase transition for some values of a^, as a function of the wall velocity v w . We chose 
AtH* = 10~ 2 , which gives f p ~ lmHz. We see that the phase transition needs not be 
too strong, i.e., with a^r > 0.3 we obtain intensities above the peak sensitivity of LISA. 
This is important, since the value of a at was found to be < 1 for several models of the 
electroweak phase transition [33]. For ~ 1 we get h 2 Q p as large as ~ 10~ 9 . 




Figure 7: The intensity of gravitational radiation from turbulence at a phase transition 
with AtH* = 10~ 2 and T* = lOOGeV, for a N = 0.03, 0.1, 0.3, 1, and 3 (from bottom to 
top). The dashed line corresponds to the peak sensitivity of LISA. 

Since only Jouguet detonations are in general considered in the literature on GWs, the 
efficiency factor is in principle overestimated with respect to real detonations. However, 
as pointed out in Ref. [22], the value of the Jouguet detonation efficiency factor kj («tv) is 
underestimated in the literature, due to a missing factor in the original paper [7J. This 
mistake compensates the fact that kj (a^) is larger than the efficiency factor n(aN,v w ) 
for weak detonations. The effect of this compensation on the intensity of GWs is shown 
in Fig. [8], where we plotted Q p as a function of v w for ct^ = 0.3, together with the correct 
(upper dashed line) and the wrong (lower dashed line) values for the Jouguet detonation 
case. The efficiency factor for weak detonations lies between these two values. Notice, 
also, that supersonic deflagrations give the largest GW amplitudes, and that even subsonic 
deflagrations can give intensities comparable to those of detonations. 

6 Conclusions 

We have studied the motion of phase transition fronts in a first-order cosmological phase 
transition, focusing on the energy injected into bulk motions of the plasma, which is a 
relevant quantity for the generation of cosmological relics such as gravitational waves. 
This issue was recently addressed in Ref. [22] for the case of spherical bubbles. However, 
the GWs are generated once the bubbles (or the shock fronts) collide, so that the spherical 
symmetry is lost. Therefore, any bubble symmetry one may assume will be just an ap- 
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Figure 8: The same as Fig. [7J for = 0.3 alone. The dashed lines indicate the results 
obtained from the correct (upper) and the wrong (lower) values of kj (o/v)- The two points 
indicate the cases v w = c s and v w = t>j et (ajv), which separate subsonic deflagrations, 
supersonic deflagrations, and detonations. 

proximation. In order to study the dependence on the wall geometry, we have considered 
bubble walls with spherical, cylindrical, and plane symmetry, for all the possible hydrody- 
namic propagation modes, namely, subsonic deflagrations, supersonic deflagrations, and 
detonations. 

We have seen that the region around the wall in which the energy spreads can be 
rather different for each wall geometry. In particular, for planar walls the region is larger, 
since the energy spreads in only one direction. For the strongest processes allowed, i.e., 
Jouguet deflagrations and Jouguet detonations, the kinetic energy of the fluid is, for the 
three geometries, concentrated in a thin region around the wall. This is because these 
processes produce a stronger disturbance of the plasma than weak processes; thus, the 
injected energy is larger and, hence, more difficult to distribute. 

The efficiency factor k, i.e., the part of the injected energy which goes into bulk 
motions (relative to the released vacuum energy), has a rather weak dependence on the 
wall geometry. This is an important result, since the walls can take arbitrary forms after 
colliding. The dependence on the wall geometry is stronger for small wall velocities, 
£ w < 0.2. For small velocities, however, the efficiency factor is small and will not play a 
relevant role in the cosmological consequences of the phase transition. Thus, it is clear 
that, for applications, it is convenient to consider planar walls, for which we have obtained 
exact analytical expressions for k (alternatively, one can use the numerical fits given in 
Ref. [22J for the spherical case). 

The efficiency factor peaks for supersonic deflagrations, which are Jouguet processes. 
Thus, k is in general sizeable for fast (close to the speed of sound) subsonic deflagrations, 
and for the slowest detonations, i.e., those solutions which are close to the Jouguet velocity. 
On the other hand, k can decrease considerably for fast detonations, and vanishes for v w — > 
0. We have applied the results for the planar case to the estimation of the gravitational 
wave signal from turbulence at the electroweak phase transition. The GW amplitude 
peaks for supersonic deflagrations. It is interesting to notice that subsonic deflagrations 
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can give efficiency factors larger than those given by detonations. Although the amplitude 
of the gravity waves may depend (according to the generation mechanism) on the wall 
velocity as well as on the efficiency factor, we have seen that subsonic deflagrations can 
produce GWs of intensity comparable to that of detonations. 
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